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We study two ions confined in a Penning trap. We show that electronically highly excited states 
exist in which an electron is delocalized among the two ions forming a giant molecule of several 
micrometer size. At energies close to the top of the Coulomb barrier these molecular states can 
be regarded as superpositions of Rydberg states of individual ions. We illuminate the possibility 
to observe coherent charge transfer between the ions. Beyond a critical principal quantum number 
the electron can coherently tunnel through the Coulomb barrier to an adjacent doubly charged ion. 
The tunneling occurs on timescales on which the dynamics of the nuclei can be considered frozen 
and radiative decay can be neglected. 



PACS numbers; 34.70.+e,32.80.Rm 

As the simplest conceivable molecule, singly ionized 
molecular hydrogen is formed out of two protons and 
a single electron. The electron is delocalized among the 
protons tying them together with an equilibrium distance 
of just a few Bohr's radii (ag). In 2000, Greene et al. 
complemented this class of conventional homonuclear bi- 
nary molecules by predicting a new type of ultra- long- 
range molecules [l| in which a Rydberg atom forms a 
molecular entity with a ground state atom. Unlike in 
the former case the electron is delocalized over two nuclei 
which are separated by 1000 or more oq. More recently, it 
has been shown that Rydberg atoms can form molecular 
states even if their electronic wave functions do not over- 
lap U . The atoms are held together by the multipole- 
interaction of their charge distributions giving rise to 
equilibrium distances of again several thousand oq. An- 
other class of 'molecular species' which however can only 
exist in a trap is constituted by trapped ions [1, Q • At 
sufficiently low temperatures the mutual Coulomb repul- 
sion which is compensated by the trapping potential leads 
to the formation of a Wigner crystal with the interparti- 
cle spacing being typically of the order of /im. The ions 
can be cooled to their motional ground state such that 
the low energy excitations are small oscillations around 
their equilibrium positions. These phonons are reminis- 
cent of vibrational molecular excitations. However, un- 
like for the previously mentioned molecular species there 
are no delocalized electrons in the system. In this work 
we show theoretically the possibility to create delocalized 
electronic states among ions which are confined in a Pen- 
ning trap d, Q . Such molecular states possess a spatial 
extension of several micrometers thereby constituting a 
novel class of trap induced quantum objects. 

We focus on a setup consisting of two ions in a Penning 
trap which is constituted by an electric quadrupole and 
a homogeneous magnetic field. In a high magnetic field, 
which we consider here, the electronic motion is essen- 



tially one-dimensional and channeled along the magnetic 
field axis. In particular, we investigate states whose en- 
ergy lies just below the Coulomb barrier separating the 
ionic charges. Such states can be thought of as being a 
superposition of states in which a single electron occupies 
a Rydberg state that can be associated with an individual 
ion. Along this line we study coherent charge transfer be- 
tween two ions, i.e. the reaction -|- + A^^ , 
which can be understood as a coherent tunneling process 
between the Coulomb wells. The transfer is initiated by 
transferring the electron of the singly charged ion to a 
high lying Rydberg state via a laser. We find tunneling 
rates of the order of hundred MHz. Over the character- 
istic time the ionic motion can be considered frozen and 
radiative decay plays no role. We discuss the dynamics 
of the charge transfer as a function of the parameters of 
the Rydberg excitation laser and the tunneling rate of 
the electron. 

The Penning trap is constituted by a homogeneous 
magnetic field B = iJe^ of strength B and an 
electric quadrupole field with the potential $(r) — 
—(3 [x^ + — 2z^] with j3 being the field gradient. In the 
symmetric gauge, the vector potential is A(r) = x r 
and the motion of a particle of charge Qe and mass M 
in this trap is governed by the Hamiltonian 



H(p,r;0,M) = ^ |p - <3eA(r)] V Qe »(r) (1) 
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Here = xpy — ypz = —ihd^ is the z-component of the 
total angular momentum, and we have introduced the 
radial ujp = y/{2/M)[{Q^e^B^)/{8M) - Qe(3] and longi- 
tudinal LOz = a/ {Qe(3)/M trap frequencies as well as the 
cyclotron frequency Uc = {QeB)/M. In order to provide 
radial confinement Up has to be real, which is achieved 
beyond the critical field strength Bi ~ ^ (8M/3) / (Qe). 
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We intend to operate the Penning trap in a regime where 
Wp > Wz, thus forming a prolate trap ^7*1. This configu- 
ration, in which the trapped ions line up in a chain is 
achieved for B > B2 = [M (3) / (Qe). 

In the following we will focus on a scenario in which 
two ions that are of the same atomic species, but carry 
different charges, are placed in the trap. The fea- 
sibility of preparing and cooling two ion systems in 
a Penning trap has been experimentally demonstrated 
in Ref. [7|. The Hamiltonian for the case in which 
the two ions are in their electronic ground state reads 
H = H(Pi, Ri; Qi, M) + iJ(P2, R2; Q2,M) + Vci\Ri - 
R-2|; Qi, (92) with the Coulomb interaction Vc(|Ri — 
R2|;Qi,Q2j = (QiQ2e2)/(47reo|Ri - R2I). We assume 
that the ionic motion has been cooled to the ground state 
such that the ions perform harmonic oscillations around 
their equilibrium positions, which are given by 



R, 



■0,i; 



k)Qi + {2-k)Q2 

(gi + Q2)2/3 



(2) 
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with the characteristic length ( = [e/(27reo/3)] 
and k — 1,2. The transverse oscillation fre- 
quencies are given by the respective values of 
the frequencies of the longitudi- 



LUp whereas 
nal modes 

Q2)) 



evaluate to 



(2e/3)/(M(gi 



Ql + 6Q1Q2 + Ql± VQf + 14QfQfTgf . 
Throughout this work we will consider ^"^Ca with 
Ql = 2 and Q2 — ^ confined to a Penning trap with the 
parameters _B = 10 T and /3 = 2 x 10^ V/m^. In this case 
we find tjpj = 2tt x 1.85 MHz, LUp^2 = 27r x 3.78 MHz, 
= 27r X 0.83 MHz and Wz,+ = 27r x 1.44 MHz. 
We furthermore have the characteristic length 
= 11.29 /im and consequently an interionic dis- 



tance of 13 = 
Bohr's radius. 



.14 /im = 153 X 10 flo with oq being 



The idea of creating delocalized electronic states in this 
setup goes as follows: The valence electron of the singly 
charged ion is excited by a laser to a high lying electronic 
(Rydberg) state. Here the electron is exposed to the po- 
tential of two doubly charged ions. The electron can hop 
as soon as its wave function overlaps with the empty Ryd- 
berg orbital of the second ion. By this electronic states 
that extend over both ionic cores are formed and the elec- 
tronic charge can be transferred. The size of a Rydberg 
orbit is approximately given by J'Ryd — (2«o''^^)/Q with 
n being the principal quantum number. Requiring r^iy^ 
to be half the equilibrium distance D between the ions 
one finds that delocalized electronic states are expected 
to occur beyond n > Uct = y/D/{2aQ) « 277 for the 
above parameters. Indeed, we will see later that this is 
the case even for much smaller values of n. 

In order to study the process more quantitatively, we 



consider the electronic Hamiltonian 

i/ei - i?(p,r;-l,m) (3) 
+Vc{\r - Roal; 2, -1) + Vc{\r ~ Ro^2|; 2, -1). 

The ions are frozen at their equilibrium position as we 
are interested in the electronic dynamics on a shorter 
timescale than the ionic motion. Moreover, the electronic 
spin is not considered as it merely gives rise to a constant 
energy offset. In Eq. ([3]) we have approximated the inter- 
action potential between the ionic cores and the electron 
by a sum of two pure Coulomb potentials of charge 2e. 
This is a simplification, since in general the true poten- 
tial will deviate from the Coulomb potential at small dis- 
tances due to the presence of the inner electronic shells. 
For Rydberg states this gives rise to the quantum defect 
of states with low angular momenta. The general fea- 
tures of the charge transfer, however, will not be affected 
by the actual value of the quantum defect and an ex- 
perimental implementation requires in any case a careful 
spectroscopic analysis. 

In the following, we will motivate an adiabatic sin- 
gle channel approximation in order to derive an effective 
one dimensional electronic Hamiltonian. To this end it 
is instructive to consider the part of -ff (p, r; — 1, m) in 
Eq. Q which governs the transversal electronic dynam- 
ics, providing harmonic confinement with the frequency 
np = v^(2/to) [(e2B2)/(8m) + e/3] = 27r x 140 GHz. 
The eigenstates of this potential are given by the two- 
dimensional harmonic oscillator states Qu,iiip,4>) with 
the energies E^,fj. — 2hilc^M + hVLp {2u + -|- 1), where 
ric = {eB)/m. We will see a posteriori that in the 
regime of interest the energy gap between the states 
Q^^p{p,4i) and Qi,±i.p,{p,<t>) is much larger than the en- 
ergy of the longitudinal electronic motion. Therefore we 
can neglect coupling between oscillator states with dif- 
ferent V. For fixed states with different p are only 
coupled if the interatomic axis does not coincide with 
the z-axis. For infinitely heavy ions in equilibrium this 
is not the case. In the case of finite mass, the corre- 
sponding couplings will lead to a modification of the 
electronic wave function depending on the positions of 
the ions which has to be taken into account by aver- 
aging over the ionic probability density. The dominant 
part of this averaging however will arise from the con- 
figuration where both ions are located on the z-axis. 
We therefore can neglect couplings between states with 
different /i and still obtain a qualitatively correct pic- 
ture. In our approximation the transversal electronic mo- 
tion is frozen in the oscillator ground state 0o,o(Pj0) — 
(\/7r/9o)"^exp [-(p//9o)^/2] with po = ^/h/{mQ.p) and 
we use the ansatz ^'(p, z, cf)) — 60,0 (p, <f>)4'{z) for the elec- 
tronic eigenstates. Such quasi one-dimensional scenario, 
in which the electronic motion is channeled along the 
magnetic field axis is typical for atoms and molecules 
(sj in strong magnetic fields. It has also been reported 
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FIG. 1: Electronic potential if the ions are frozen at their equi- 
librium positions. The electronic potential supports bound 
states which are delocalized between the two nuclei. For too 
large excitation energies the saddle point due to the inverted 
harmonic potential leads to field ionization. 

for electron-hole pairs in a semiconductor at high mag- 
netic fields [9] and for long-range molecules formed by 
a ground state and a Rydberg atom The effective 

Hamiltonian for the longitudinal dynamics is obtained 
H\\ ^ J d(f>J dpp&^ Q{p,(l))Hcieo,o{p,(j)). Neglecting 
a constant energy offset one obtains 

xerfcf^^^V-™^^'^' (4) 
V Po / 2 

where D,^ = 2^(e/3)/m = 27r x 189 MHz is the 'fre- 
quency' of the inverted parabolic potential and erfc(a;) 
is the complementary error function. The correspond- 
ing longitudinal potential is presented in Fig. [TJ If the 
electron is far away from either ionic core the potential 
is approximately a sum of two one-dimensional Coulomb 
potentials Vid = —e'^ / {2TTeo\z\) centered at the respective 
ionic equilibrium positions. For small excitation energies 
the potential wells are well separated. If the excitation 
energy surpasses Sdoi « -(2e2)/(7ren£') the wells merge 
and hence electronic states which extend over both ionic 
cores are possible. Expressing E^ci in terms of the prin- 
cipal quantum number n one finds this merging to oc- 
cur approximately at ridoi = (e/fi,)-^/ (mD) / (167reo). For 
the employed parameter set one obtains ridoi ~ 196.14 
which is significantly smaller than the principal quan- 
tum number we obtained earlier by invoking the argu- 
ment of overlapping wave functions. Further inspection 
of the potential reveals that the inverted parabolic po- 
tential due to the Penning trap eventually leads to field 
ionization which is shown in Fig. [1] The position of the 
saddle (zgd > ^0,2) is given by Zsd = 7^0,2 where 7 > 1 
is the solution of the equation (7 — 1)~^ + (7 + l/2)~^ — 

(77-0 and a = [{\/2)mnlZl^] [eV(4^eo|^o,2|)] = 
(l/2)Qf/(Qi + Q2Y = 4/9. This equation is solved by 
7 = 2.12 and thus the saddle point is for our param- 



FIG. 2: Sketch of the energy level structure close to ridoi. 
Here the inverted parabolic potential can be neglected and 
the states can be characterized by their symmetry property 
under reflection at Zsym. The energies of these pairs of ger- 
ade/ungerade states are indicated by dashed/solid lines. The 
pairs are almost degenerate with a small energy splitting m. 
Wave functions which are localized in one of the two wells 
{\Ei/r)) can be created by linear superpositions of the gerade 
and ungerade states. Two such wave functions are exemplar- 
ily shown. 



eters located at Zsd — 217523 aq. From this we find 
that classically ionization is expected to occur for prin- 
cipal quantum numbers beyond nion = 1.456y'.Zo,2/ao = 
(0.21e2/3^)/(&^/^/3i/6) 466. Hence there is an en- 
ergy window ridoi < n < Tijon in which delocalized elec- 
tronic states exist. Here the escape of the electron from 
the trap is prevented by its attraction to the doubly 
charged ion. For the given parameters the width of this 
window is about 60 GHz. 

We now focus on the electronic states which lie close 
to the top of the Coulomb barrier. In this regime the 
influence of the anti-confining potential can be approxi- 
mately neglected and the problem assumes a particularly 
symmetric appearance. We obtain the eigenstates of in- 
terest by numerically integrating the Schrodinger equa- 
tion of the Hamiltonian (|4]) with = in the vicinity of 
ndc\- With respect to reflections at the symmetry point 
^^sym — {Zo,i + -2^0,2 )/2 we can define gerade and unger- 
ade states. In Fig. [2] we show a sketch of the energy 
level structure in the vicinity of ridci. Gerade {\Eg)) and 
ungerade states {\Eu)) are energetically split by the en- 
ergy 2k = Eu — Eg with K being the tunnel coupling. The 
splitting increases with increasing degree of excitation. 
States which are localized in the left {\Ei)) or right well 
{\Er)) can be created by the linear combinations \Ei) = 
il/V2)[\Eg) + \Eu)] and \Er) = (1/^2) [ji^^) - |£;„)], 
respectively. They are not stationary but evolve under 
the Hamiltonian i?tunnGi = 1^'/) {Er\ + h.c. which de- 
scribes Rabi oscillations of the electronic density between 
the left and right well at a rate k. As an example we pro- 
vide two sets of gerade/ungerade eigenstates which were 
obtained for the above parameters: Their energies corre- 
spond to the principal quantum numbers = 195.6848, 
ni = 195.6917 and = 196.0161, nl = 196.0662 
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FIG. 3: Cut through the potential curves along the normal 
coordinate Za with Zh = 0. The equilibrium positions of the 
initial configuration (A"'"^ + ^4"*") do not coincide with those of 
the final configuration (A^ + A^^'). As a consequence there 
is a relative detuning A between the initial and final state 
which has to be accounted for during the laser excitation. 



(E(n) = -4 X 13.6eV7i"2)^ ^hc fact that these values 
are smaller than yijci indicates that the electron in the re- 
spective states indeed tunnels through the barrier which 
separates the two Coulomb wells. The tunnel coupling 
in the two cases evaluates to Ki/h= 2tt x 12.3 MHz and 
K2/h = 27r X 87.5 MHz. Since these rates are significantly 
larger than the oscillation frequencies of the ionic cores 
it is indeed justified to use the picture of a frozen ionic 
motion. Moreover, the rates are much larger than the 
radiative decay rate of Rydberg states which is usually 
smaller than 1 MHz. Therefore the described process can 
be considered coherent. 

We are now in position to study the charge transfer 
process yl"*""*" + — > + A"*""*". To this end we employ 
two laser fields coupling the ground state of the right ion 
\Gr) resonantly to the state \Er) with a Rabi frequency 
ri2- Due to the tunnel coupling the state \Er) evolves 
into the state \Ei) which is then resonantly coupled to 
the ground state |G/) of the left ion with Rabi frequency 
fli. While initially the two ions resided in their equi- 
librium position this is no longer true for the final state 
as here the charges are interchanged. The new equilib- 
rium positions can be calculated by Eqs. However, 
if the ions are frozen, their energy in the final state, i.e. 
j4+ -|- A"*""*", is higher than in equilibrium. Hence, one 
has to account for this energy off-set A in order to make 
the transition \Ei) \Gi) resonant. This is shown in 
Fig. [3] where the potentials of the ground state configu- 
rations A++ + A+ and A+ + A++ arc depicted. The in- 
tersections shown are made along the normal coordinate 
Za = a-Zi + olj^Zi at Z\, — ~aj^Z\ -I- oi-Zi set to zero. 

The coefficients are given by a± = \J (1/2) ± 3/(2-\/7^. 
The energy gap between the two configurations evaluates 
to A = (e/3^^)/(2 X 3^/'^) which yields for our parameters 
A = 27r X 844 GHz. 

With the two lasers turned on and choosing = ^2 — 
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FIG. 4: Time dependence of the probability to find the elec- 
tron with the left ion for 5 — 0.2, 0.5, 2.5 and the asymptotic 
curve for 5^1. Initially the electron is with the right ion. 
The time-evolution is governed by the Hamiltonian ([5]). 

f2 the electronic dynamics is governed by the Hamiltonian 

iftunncl = n \Gl) {Ei\+n \Gr) {Er\ + K \Ei) {Er\ + h.C.,(5) 

where \Gr/i) is the ground state of a singly charged ion 
at the left or right position in the trap. The probability 
of finding the electron in the left potential well is given 
by 

Picft(i) = ^ _^^g2 [S'^ {cos{e -Kt) - cos{e+Kt) f 

-f (e+ sin(e_Kt) — e_ sin(e_|_Kt))'^] (6) 

with e± = y'(l + 2(5± Vl + 4(52)/2 and S = n/n. In 
Fig. 0] this probability is shown for three different val- 
ues of S. For S = 0.2 the electron mainly stays with the 
right ion over the time interval shown. As 5 increases 
the likelihood to find the electron in the left well in- 
creases significantly for small times. In the limit 6^1 
one finds picft(t) = (1/2) (1 — cos(Kt)) and the electron is 
found with certainty to be with the left ion at the times 
t„i = (2m + 1)tt/ K. Once the transfer is performed and 
the lasers are switched off the system is in a highly excited 
state as can be seen from the potential surfaces presented 
in Fig. [3] A successful charge transfer will be accompa- 
nied by strong oscillations of the ions. After cooling the 
ionic motion the interchange of the charge can be mon- 
itored by fluorescence imaging on a frequency which is 
only resonant on a transition of the singly charged ion. 
The above treatment is only valid for time intervals over 
which the ions don't move. Moving ions would turn the 
tunneling rate into a dynamic quantity. 

In principle the presented scenario can be extended to 
more than two ions. An electron could then form states 
which are delocalized over an entire ion chain. Such 
system would constitute an implementation of the tight- 
binding model. Future directions of this work are to cre- 
ate and to study one- and two-dimensional realizations 
of such system and to include more than one electron in 
the dynamics which would give rise to Hubbard-models. 
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